ON THE EVOLUTION OF THE LIGHT ELEMENTS 

1. D, ^He, AND ^He 
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; ABSTRACT 

^ ■ The light elements D, ^He, ^He, and ^Li are produced in big bang nucleosynthesis and 

O ' undergo changes in their abundances due to galactic processing. Since one may observe 

^ ! most of these elements only in contemporary environments, knowledge of the interven- 

V/-^ I ing evolution is necessary for determining the observational constraints on primordial 

Q^ ■ nucleosynthesis. Chemical and stellar evolution model dependences in light element 
evolution are systematically investigated via a comparison of 1460 possible chemical 

^! evolution scenarios and of stellar nucleosynthesis yields, all of which have been selected 

O I to fit solar neighborhood C, N, O, and Fe abundances as well as the observed local 

• gas density and gas mass fraction. The light element evolution and solar system yields 

ci ■ in these models are found to span a wide range, explicitly demonstrating the model 
dependence. The range of model dependence for D, ^He, and ''He solar abundances is 



calculated, and its sensitivity to the heavy element constraints is noted. The chemi- 



;_i ■ cal evolution contribution to the uncertainty in the observed primordial light element 

^ abundances is estimated, and the effects of this uncertainty on big bang nucleosynthesis 

results are discussed. The predictions for the light elements are found to be correlated; 
the extent and physical origin of these correlations is discussed. D and ^He evolution is 
found to have significant model dependence, however, the dominant factor determining 
their solar and interstellar abundances is their primordial abundance. In addition, ^He 
is found to be very sensitive to the details of processing in low mass stars. ^He yields 
are shown to be very model dependent; in particular, both the introduction of mass 
loss and the possibly very high ^He stellar yields in the poorly understood mass range 
of ~ 8 — 12 Mq can lead to large enhancement of ^He production and can lead to large 
slopes of AY/ AN and AY/ AO. It is found that the inclusion of secondary nitrogen 
leads to only a small distortion in the low metallicity Y— N relation if there is also a 
significant contribution of primary N, as required by observations. 

Subject headings: Galaxies: evolution — nuclear reactions, nucelosynthesis, abundances 
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1. Introduction 

A central element of the hot big bang cosmo- 
logical framework is the prediction of light el- 
ement (D, ^He, ^He, and '^Li) synthesis during 
the epoch of big bang nucleosynthesis (hereafter 
BBN). This prediction is tested against observa- 
tions of the light elements abundances, and the 
agreement between theory and observation pro- 
vides major empirical support for the big bang. 
The comparison between prediction and mea- 
surement is not immediate, however. Many of 
the observations are of environments that are 
contemporary to the present, or nearly so, and 
galactic astrophysical production and destruc- 
tion of the light elements may have changed their 
abundances. An understanding of these galactic 
processes is thus needed to in order to accurately 
interpret the light element observations and so to 
realistically constrain BBN. The study of the nu- 
clear history of galaxies is the domain of galactic 
chemical evolution; here we will re-examine the 
galactic chemical evolution of the light elements. 

Models for galactic chemical evolution were 
first constructed in the 1970's and for the most 
part studied the nucleosynthesis history of the 
heavy elements. Since such pioneering work as 
Cameron & Truran ( |1971| ), Talbot k Arnett 
( |T97ll ), and Tinsley ( |l980| , |197^ ), chemical evo- 
lution has come to do a surprisingly good job 
of reproducing many solar abundance patterns 
and even many abundance histories as observed 
in and around the solar neighborhood (see, e.g., 
the recent work of Timmes, Woosley, &; Weaver 
1995; Matteucci & Francois 1989| ; and Ferrini 
et al. 1992 ). These models also give reason- 



able predictions for stellar luminosity evolution 



(see, e.g. Tinlsey |1980|; Colin Sz Schramm |1993[) , 
and underpin the study of nucleocosmochronol- 
ogy (e.g. Thielemann et al. 1993). 

The chemical evolution of the light elements 
has been studied since the work of Cameron & 
Truran ( |1971 ), and has been used to determine 
primordial abundances since the work of Reeves, 



Audouze, Fowler, & Schramm (1973) and Au- 
douze & Tinsley ( |1974D . These seminal papers 
made a critical comparison of light element abun- 
dances from deuterium to boron with the vari- 
ous cosmological and astrophysical processes sug- 
gested for their production. Since this early work 
there have been many calculations of light ele- 
ment chemical evolution, usually focusing on a 
few nuclides (see references in §§ ^.3j - 2".5.1.| ). In 
this paper we will again take a comprehensive 
look at all of the light elements to see how con- 
sistent a picture of their joint chemical evolution 
can be made. 

Chemical evolution models are physically based 
on important mass conservation principles, and 
have been refined considerably over the years. 
Unfortunately, though, chemical evolution is still 
a somewhat ambiguous enterprise. One should 
always be mindful of the uncertainties in the 
conventional prescription for these calculations. 
These uncertainties stem from unknown charac- 
ter of star formation rates and mass distribu- 
tions, galactic inflow and outflow of gas, and from 
unresolved model dependent effects in stellar evo- 
lution nucleosynthesis calculations. Lacking a 
first principles calculation of these phenomena, 
we are forced to model them in a phenomenolog- 
ical or schematic way. 

As a result of these uncertainties, the com- 
monly adopted framework for chemical evolution 
contains features that are not well understood 
in a theoretical and/or empirical way. Thus 
within this framework, there are a great many 
ways to construct specific models, and work has 
been done using a diverse variety of the possible 
scenarios. Encouragingly, the overall framework 
has been quite successful despite these difficul- 
ties. And indeed chemical evolution is a always 
a part of comparison of nucleosynthesis calcula- 
tions with data (though sometimes its role is only 
implicit). 

The effects of chemical evolution model de- 
pendences have been investigated to varying de- 
grees. However, the larger studies of these effects 
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(e.g. that of Tosi 198^ ) did not concentrate on the 
hght elements, and work on hght element evo- 
lution has not examined model uncertainties as 
broadly. There has not yet been an emphasis on 
the full range of uncertainty in a comprehensive 
description of light element chemical evolution. 

In this paper we construct a chemical evolu- 
tion framework for the evolution of all of Ithe 
light elements, and we obtain estimates of the 
range of model uncertainties. We will follow the 
history of all of the light elements, D through 
^^B. We also include the non-cosmological "heavy'' 
elements ^^C, ^^N, ^^O, and ^^Fe, which serve 
as tests of potential models. We systemati- 
cally examine a wide variety of model assump- 
tions, examining many chemical evolution pre- 
scriptions and parameterizations (in the spirit of 



Tosi 1988 ), as well as including many suggested 
galactic production and destruction mechanisms 
for each of the light elements (in the spirit of 



Reeves, Audouze, Fowler, & Schramm 1973 ). 

Our strategy is to choose a broad chemical 
evolution framework, allowing for many features 
that authors have considered in examining the 
histories of specific light elements. These fea- 
tures are chosen to be inclusive, without regard 
for their compatibility. With these we form a 
very large (~ 1200) set of models which include 
all combinations of these features. We then con- 
strain these models to fit observed solar neigh- 
borhood characteristics, apart from those of the 
light elements. This leaves ~ 120 potential mod- 
els which are viable descriptions of solar neigh- 
borhood chemical evolution (again, light element 
behavior aside) . To each of these successful mod- 
els we add a range of 18 variants of light element 
evolution. We then examine the predictions of 
this large battery of light element evolutionary 
schemes. Specifically: (1) we see which mod- 
els fit light element observations; (2) we estimate 
the chemical evolution uncertainty and discuss 
its influence on the use of BBN observational 
constraints; and finally (3) we note correlations 
between the evolution of different light element 



nuclides. As part of this analysis we will consider 
as well the uncertainties in the solar system con- 
straints themselves. 

In addition to these considerations of the whole 
set of the light elements, we will also address spe- 
cific issues particular to each of them. For exam- 
ple, recent observations of high-redshift quasar 
absorption line systems have approached the pre- 
cision needed to give meaningful D abundances. 
This is an exciting prospect indeed, as these high- 
z environments are much more pristine than con- 
temporary ones, and so give a new window to a 
much earlier epoch in D evolution (though still 
> 1 Gyr after BBN). The initial reported abun- 
dances by Songalia, Cowie, Hogan, & Rugers 



(1994) and by Carswell, Ranch, Weymann, Cooke, 
& Webb ( |1994 ) were tantalizingly high if taken at 
face value (see, e.g. Kernan & Krau ss |1994| ; Casse 
& Vangioni-Flam 1994; Steigman |1994| ). How- 
ever, the observers themselves have been care- 
ful to caution that the absorption line systems 
may be subject to contamination, and so the de- 
rived D abundances are best understood as up- 
per limits to the actual levels of D/H, a caution 
strengthed by recent reports of a low D abun- 
dance in a different line of sight (Tytler &; Fan 
1994). Clearly, the present situation is in flux and 



more observations are needed to determine the D 
abundance in these systems accurately; but just 
as clearly, this technique is quite promising and 
holds the key for a very accurate D abundances 
at very early epochs. We will examine D evo- 
lution in detail in this paper, and comment on 
the role of it chemical evolution in interpreting 
the quasar absorption line system observational 
program. 

These new constraints on D evolution also call 
attention to conventional accounts of the evolu- 
tion of ^He. Recent improved observations of ^He 
in Galactic H II regions can be puzzling in light of 
preceding discussions of this element. The data 
suggest that ^He shows spatial variations, and 
that the abundance of ^He today in the ISM is 
not significantly larger than that at the forma- 



3 



tion of the solar system — indeed, in some places, 
^He is even seen in (possibly) smaller amounts. 



This may call into question (Hogan 1995) the 
conventional wisdom about ^He production in 
low to intermediate mass stars (Iben & Truran 
1978| ). On the other hand, the high ^He abun- 



dance seen in planetary nebulae (Rood, Bania, 
& Wilson \m% ) seems to support the notion of 
production in low mass stars at a level roughly 
consistent with that predicted by Iben & Truran 



1978| . Thus the H II region data may more likely 
point to a need to re-examine ^He evolution in 
high mass stars (Olive, Rood, Schramm, Truran, 



& Vangioni-Flam 1994). At any rate, the issue 
of low mass stellar production of ^He is crucial 
and we will address it explicitly. 

The much more abundant helium isotope, ^He, 
is important because it is so powerful. No- 
tably, an accurate assay of the ^He content of 
the universe provides a bound on N,y, the num- 
ber of light neutrino families, as first shown by 
Steigman, Schramm, and Gunn ( 1977| ). Thus 
there are continual improvements in the accuracy 
of theoretical predictions for ^He (e.g. Dicus et 
al.p82| ; Walker et al. p9T| ; Kern an |T993t Seckel 
19941; Fields, Dodelson, & Turner |1993|) , as weh 
as in the observed abundances (Skillman et al. 
1994| ). A pressing issue now is the systematic 
uncertainty in these measurements (e.g. Copi, 
Schramm, &: Turner 1995; Sasselov & Goldwirth 



1995| ). There is also discussion of the correct 
means of extrapolating the primordial '^He abun- 
dance, an issue in which chemical evolution plays 
a crucial role. 

In this paper we will describe our framework 
for constructing models of light element chemi- 
cal evolution, and summarize our overall results 
on model dependences. We will then discuss in 
more detail our results for D, ^He, and ^He evolu- 
tion and its uncertainties. In a forthcoming work 



(Fields 1995, hereafter Paper II) we will discuss 
our results for Li, Be, and B evolution, with par- 
ticular focus on the range of available models for 
cosmic ray nucleosynthesis of these elements. 



2. Light Element Evolution 
2.1. BBN Predictions 

The BBN predictions are crucial to our con- 
siderations, as the primordial yields set the initial 
conditions for the chemical evolution calculation. 
The basic BBN calculation itself has not changed 
dramatically since the calculation of Wagoner, 
Fowler, & Hoyle (|1967|) . For reviews we refer 



the reader to, e.g., Schramm & Wagoner ( 1977 ), 
Boesgaard & Steigman (I1985D , and to Smith, 
Kawano, &: Malaney ( 1993|). For recent calcula- 
tions see Walker et al. (|l99ll) , Krauss & Kernan 
(1994), and Copi, Schramm & Turner ( |1995 ). 

We consider only the simplest and "standard" 
case, namely that of BBN in a universe that 
is spatially homogeneous in all particle species. 
Some models with baryon inhomogeneities are 
still allowed by the data, but even these more 
complicated cases lead to the same basic conclu- 
sions, as discussed in, e.g., Malaney &: Mathews 
([I993D and Thomas et al. (|1993D . 

In the standard calculation, the only free pa- 
rameter is r] = ns/n^, the baryon-to-photon ra- 
tio. In Figure |^ we plot the light element abun- 
dances as a function of 77; for reference, Copi, 
Schramm, & Turner ( |1995| ) found concordance 
between theory and data for the range 



2.4 X 10"^° < 7? < 6 X 10 



-10 



(1) 



which we have indicated on the plot. 

For our purposes there are several salient fea- 
tures to notice. First, the ^He abundance does 
not depend strongly on r/, with variations only 
at the 10% level; however, we will be interested 
in the ^He abundance to this level. Even more 
evident is the decline in D and '^He abundance 
with increasing ry, allowing a large range of ini- 
tial values, particularly for D which varies more 
strongly with 77. This point is important, as we 
will find that the a requirement of significant de- 
pletion of initially high D down to the solar and 
interstellar abundance is a serious constraint on 
chemical evolution models. 
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As we will see, the combined abundance D+^He 
is a most useful constraint on rj. Note that for 
the concordant range of rj that this sum is dom- 
inated by the contribution from D. Knowledge 
of the primordial level of D is tantamount to 
knowledge of D+^He; thus the D abundances ob- 
tained from quasar absorption line systems pro- 
vides a constraint on this sum. Thus an accu- 
rate quasar measurement would provide a crucial 
complement to the solar system and interstellar 
D and '^He abundances. 

2.2. Chemical Evolution Framework 

The basic theoretical framework and observa- 
tional constraints on models of chemical evolu- 
tion is reviewed in Tinsley ( |198C| ). We will, for 
the most part, adopt her notation. Our goal is to 
trace the history of the Galaxy's isotopic content. 
Ideally, we wish to model formation of stars in 
inhomogeneously distributed gas clouds, their re- 
turn of material during their lifetime to the ISM 
via winds, and finally their deaths, which occur 
in stochastic bursts and return processed and un- 
processed material to the ISM. In principle, all 
of these processes can depend on many variables, 
e.g. the metallicity of the gas, the temperature of 
the gas, and the strength and configuration of the 
galactic magnetic fields, to name a few. 

However, this depiction is not computation- 
ally feasible. To reduce to problem to a tractable 
one, we will employ several simplifying assump- 
tions: (1) the Galaxy (or at least the solar neigh- 
borhood) is spatially homogeneous; (2) stars are 
born according to a given mass distribution, and 
their lifetimes are a function of mass only; and 
(3) all of the stellar ejecta are returned at once 
at the time of the star's death, and instanta- 
neously mixed in the interstellar medium (ISM). 
These assumptions, though strong, are adequate 
for many applications, and in any case are nec- 
essary to proceed. 

We will employ the standard convention for 
the stellar creation rate C, which gives the num- 
ber of stars born (per pc^) per Gyr, in dividing 



it into a star formation rate (hereafter SFR, ^) 
and an initial mass function (hereafter IMF, (p) 



C{m, t) dm dt = 'ip{t)(f){m) dm dt 



(2) 



The SFR has units of [MQGyr^-'^pc^-'^], and gives 
the mass of ISM material going into stars per unit 
time. The IMF gives the mass distribution with 
which these stars are formed. Questions of the 
time constancy of IMF have been to answer theo- 



retically or observationally (see, e.g. Scalo |l98q) , 
although recent work (Beers, Preston, & Shect- 
19921) supports the conventional assumption 



man 



of its constancy, which we will follow. 

The literature is fraught with unfortunate vari- 
ations in the convention for IMF normalization; 
we follow Tinsley ( 1980| ) in choosing 



dm m (j){m) = 1 



(3) 



miow 



The lower and upper mass limits, miow and rriup, 
are not well constrained by data, and so will be 
parameters we will have to choose. 

We will use a one zone model for the solar 
neighborhood, treating it as homogeneous and 
instantaneously mixed. We will track the total 
disk mass in terms of its surface mass density 
fjtot, which evolves by 



ot 



dt 



f 



(4) 



where / is the rate for the infall of new material 
due to, e.g. the collapse of the disk or to accre- 
tion of extragalactic material, in units [M^gGyr"^ 
pc-2]. 

We put the basic constituents of matter as 
stars (including their remnants) and gas. 



0"tot = (Ts + fJg . 



(5) 



Of these, we will explicitly track the gas, which 
evolves according to 



do's- , ^ „ 



(6) 
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where 



E{t) 



niup 



dm nf^ (j)(m) ip {t — r^) (7) 



■nit 



is the rate at which dying stars eject processed 
material back into the ISM. The lower mass hmit, 
mt, is the stellar mass 771 for which. = t, i.e. the 
smallest mass star which dies just at time t after 
being created at t = 0. 

We will define the mass fraction 



X,, = 



C^tot 



X,, = 1 



and we write 



dt 



-ijXi + E, + fX, 



inf 



(8) 



(9) 



where 
Ei{t) 



dm mf(l){m)ip {t — r^) (10) 



is the rate of the ejection of isotope i to the ISM. 

Note that the presence of the "retarded time" 
t — Tm makes the equations unwieldy. In practice 
one often makes ip a function of Ug and/or atot, 
thus rendering the basic equations as integro- 
differential and hence without an analytic solu- 
tion We have therefore numerically implemented 
the preceding equations, taking into account their 
full integro-differential nature. Our results are all 
based on these calculations. 

For comparison with calculations using the in- 
stantaneous recycling approximation (IRA), it is 
useful to define the return fraction 



R 



dm m^^(j){m) 



(11) 



the fraction of mass going into a generation of 
stars that will eventually be ejected back into the 
ISM. Note that i? < 1 and that it is very sensi- 
tive to the choice of IMF form and mass limits. 
We will also use fi = cTg/atot for the gas mass 
fraction. 

Of the model features we have described, those 
most influential on the abundance predictions are 



the IMF, the SFR, infall, and the adopted nucle- 
osynthesis yields. Of these, the IMF, SFR, and 
infall (or outflow) are not well understood the- 
oretically, and so one must select a phenomeno- 
logical prescription for each. The stellar nucle- 
osynthesis inputs come from model calculations 
for different stellar masses (and sometimes differ- 
ent metallicities) . For the most part only a lim- 
ited number of metallicities are used, thus forcing 
chemical evolution models to interpolate and ex- 
trapolate. A better treatment would iteratively 
run both stellar models and chemical evolution 
models; such a program has recently been done 



by Timmes, Woosley, Sz Weaver (1995), who em- 
phasized heavy element yields which they calcu- 
lated self-consistently using a supernova code. 

Of the other input parameters, two are cos- 
mological. The first is of course the BBN result 
of the initial abundances; our approach to this 
is discussed below (§^). The other cosmological 
parameter is the age of Galaxy including the col- 
lapse of the disk; we will take this 15 Gyr, and 
our results are not very sensitive to this choice. 
We will take the age of the earth to be 4.6 Gyr. 
Finally, we will use the Tm vs. m relation of Scalo 



(p86|) . 

Having selected a set of models represent- 
ing the range of choices of the major model 
features — IMF, SFR, and nucleosynthesis — we 
will want to constrain this set with observational 
data. Specifically, we will test these against so- 
lar system abundances of CNOFe, as well as 
the gas mass fraction. The confrontation with 
observations will select a set of viable models 
which we will then use to test different models 
of light element abundances. Thus we will want 
to include in our initial suite of models all vari- 
ations of the major features — the IMF, SFR, 
and nucleosynthesis — which have been proposed 
to be well-suited to fit individual light elements. 
Therefore, we now will review the chemical evo- 
lution of the light elements, to identify features 
to include in the initial set of chemical evolution 
features. 
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2.3. Deuterium 

2.3.1. D Data 

Traditionally, D has been measured in the so- 
lar system and in the ISM. Important evidence 
for the solar system D abundance came as a re- 
sult of the Apollo 11 mission (Geiss &: Reeves 
1972| ), giving a measurement of the D+^He con- 



tent in the solar wind. These data are comple- 
mented by measurements of meteorites (Black 
1972| ), which allows determination of the preso- 



lar D level. The best numbers now are that the 
presolar D abundance is (D/H)0 = 2.6 ± 1.0 
(Geiss [1993| ). 

The D abundance in the ISM was first mea- 
sured by Rogerson & York ( p^73|) . It has most 
recently been observed by the Hubble Space Tele- 



scope, with Linsky et al. (|1993|) reporting (D/II)ism 



l-Slg^g- This single, remarkably precise mea- 
surement will prove to set strong constraints on 
chemical evolution models. 

Most D measurements to date have been of 
local and relatively recent cosmological epochs. 
However, exciting new observations have reported 
possible detection of D in quasar absorption line 
systems (QSOALS). These systems considered lie 
at high redshift {z > 3), and so are very young 
(1-3 Gyr). Thus one expects D to have suf- 
fered much less depletion. Indeed, initial reports 



(Songalia, Cowie, Hogan, & Rugers 1994 ; Car 



swell. Ranch, Weymann, Cooke, & Webb |1994D 
suggested the possibility of a very high D, at lev- 
els of D/H < 2.5 X 10~^. Both of these groups 
have been careful to caution that this result can 
only be regarded as an upper bound on D, as 
low column density in intervening clouds can con- 
tribute significantly to the putative D line. In- 
deed, this may be the case for these measure- 
ments, as Tytler & Fan ( |1994| ) have recently ob- 
served a different line of sight and in a high- 
redshift absorption system claim a preliminary D 
abundance of D /H ~ 2 x 10~^ . The situation may 
be more muddled still, as it has recently been 
suggested (Levshakov & Takahara 1995| ) that im- 



proper modeling of the turbulent characteristics 
of these systems could lead to very large erros in 
the derived D abundance. We will thus take the 
prudent route of not using the few D abundances 
obtained in this manner, and instead concentrate 
on better-determined the solar and interstellar 
values. 

2.3.2. D Chemical Evolution 

Deuterium was originally thought not to be 
primordial, but to be produced in the T-Tauri 
phase of stellar evolution (Fowler, Greenstein, &; 
Hoyle 1962| ). However, this process was shown 
to fail (Ryter, Reeves, Gradstajn, & Audouze 



1970 ), and Reeves, Audouze, Fowler, &; Schramm 
( 1973 ) argued that D was likely of cosmological 
origin. This argument was cemented by the work 
of Epstein, Lattimer, & Schramm ( 1976| ), who 
showed that not only does stellar burning destroy 
D, but any high energy astrophysical process de- 
stroys it as well — except the big bang. 

Indeed, as recently emphasized by Copi, Schramm, 
& Turner ( 1995| ), D is best baryometer of BBN 
in that the constraints on ry which arise from D 
observations are particularly reliable. This reli- 
ability stems from the uniquely straightforward 
evolution of D in galaxies: all of the D in the 
universe was produced in BBN, and all subse- 
quent processes destroy it. This fundamental 
fact greatly simplifies the chemical evolution of 
D, which has no sources and so is absent from 
any processed material that is returned to the 
ISM. The D abundance thus provides a clean in- 
dication of the amount of material that has not 
been processed in stars. (The galactic evolution 
of D ^He, and ^He is summarized in table |l]). 

To get a feel for D evolution, we note that in 
the instantaneous recycling approximation, the 
D evolution can be solved exactly as 



(12) 



Note that this result does not depend on the form 
of the star formation rate ip. Furthermore, since 
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the the gas fraction fi is constrained by observa- 
tion, the entire character of the depletion comes 
from R, the return fraction. This is in turn com- 
pletely determined by the initial mass function. 
Physically, the D abundance is a measure of the 
fraction of gas which remains unprocessed. Thus 
as material is cycled into stars, the D depletion 
will depend on how much (D-free) material is re- 
turned to the ISM to dilute the unprocessed com- 
ponent. 

The main chemical evolution feature to con- 
sider in order to capture the range of D evolution 
is that which determines the return fraction R, 
namely the IMF. We will want to consider the 
gamut of possible IMF forms and mass limits. 
Indeed, we will expect all the light element evo- 
lution to depend strongly on the adopted IMF. 

2.4. Helium 3 

2.4.1. '^He Data 

The solar system information on both ^He and 
D come together and are derived from measure- 
ments of the solar wind and meteorites (Geiss 
[1993| ; Black 1972); the presolar abundance is 
(^He/H)© = 1.5 ± 0.3. In the ISM, ^He is ob- 
served in Galactic H II regions, and measured via 
its hyperfine line. Although these measurements 
are difficult, they apparently show real varia- 
tions, with (3He/H)isM ~ (1 - 5) x 10"^ (Balser, 



et al. 1994| ). ^He has also been detected by Rood, 
Bania, & Wilson ( |l99^ ) in a planetary nebula, 
and is found at very high levels: ^He/H ~ 10"^. 

2. A. 2. ^He Chemical Evolution 



As noted in S2.4 



BBN yields a fairly low '^He 
abundance, but not much lower than contempo- 
rary observations; i.e. BBN theory does not leave 
a lot of room for the increase of ^He. On the other 
hand, stellar models firmly show that all initial 
D in a star is processed to ^He in the pre-main 
sequence phase when the star is fully convective. 
Some of this ^He survives to be ejected at the 
star's death, and indeed in low mass stars are 



likely to be ^He sources. Thus while D suffers 
astration and '^He is produced over time, the sum 
of the two, D+^He, is more stabile; furthermore, 
D+^He is first dominated by D and then ^He 
increases until roughly equal measures of each 
comprise the present ISM abundance. 

Yang et al. ( |1984D get a bound on (D+3He)p 
by noting that in all stars, any initial D goes 
to ^He in the pre-main sequence phase. Some 
fraction g^{m) < 1 of that ^He survives in the 
ejecta at the star's death. If we consider one 
stellar generation, and let ((73) be the average of 
53 (m) over the mass function, it may be shown 
that 



'D + 3He' 



< 



1 



(53 



H 



+ 



(Yang et al |1984D . The inequality (eq. |l^) comes 
from disregarding ^He production, and so guar- 
anteeing that X23 decreases with time. 

We may get another perspective on this point 
by seeing how the D+^He argument plays out not 
for a single stellar generation, but in the approx- 
imation of instantaneous recycling (see preced- 
ing section) for a closed box chemical evolution 
model. We have 



^X^ + R{g^)^[ Xy,+X^] (14) 



d (q-gX3) 

dt ' ' "^•"-^'•^ V2 

and when substituting from eq. |^ we have 



(IX2 + X3) 

(1^2 + X3)p 



X2 



i-(5'3 



(15) 



and so in term of observable number ratios, 

^-{93) / V \ (93) 



(D + 3He) 





^-{93) 



(16) 



D+^He increases (decreases) with time if (53) > 
1 ((53) < !)• Note, however, that because low 
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mass stars are crucial here, that the IRA is in- 
appropriate. At any rate, there is a close rela- 
tionship between D and ^He, and we see that it 
is crucial to know whether ^He is created or de- 
stroyed with time. 

The physical explanation for the different fate 
of ^He in high and low mass stars is that for 
m < 2Mq , p — p is the dominant form of hydro- 
gen burning, and the p{p,i'e^)D{p,j)^ILe chain 
makes '^He. However, for m > 2Mq, the CNO 
cycle dominates the hydrogen burning, and ^He 
is destroyed. Indeed, Iben & Truran's ( 197^ ) 
calculation for low mass stars suggests that ^He 
is produced copiously. Dearborn, Schramm, & 
Steigman ( |1986| ) find that ^He is destroyed in 
high mass stars, particularly at low metallicity. 



Woosley & Weaver ( [L994| ) give more detailed but 
qualitatively similar results. 

These conventional notions about ^He evolu- 
tion have recently been called into question as the 
interstellar ^He abundances have become more 
accurate. Until recently, the production of ^He 
in low mass stars was not considered in chem- 
ical evolution studies, e.g. that of Steigman &: 
Tosi ([19921) . Vangioni-Flam, Olive, and Prant- 
zos ( tl994| ) let the contribution to (53) from low 
mass stars be a free parameter, and used the 
Dearborn et al. ( 1986| ) high mass yields. They 
found that to fit the data require net destruc- 
tion by low mass stars, a mechanism for which 



is suggested by Hogan (|1995D . Indeed, with the 
Iben k Truran ( [l978| ) and Dearborn et al. (|l98|) 
results, Galactic ^He which starts with the es- 
poused BBN value is apparently overproduced 
in many chemical evolution models (Olive et al 



1994| ; Galli et al. |1994|) . However, in the face of 
this apparent need for '^He destruction is the ob- 
served prescence of high '^He in planetary nebulae 
(Rood, Bania, & Wilson |1992| ) at a level in rough 
agreement with Iben & Truran's calculation. 

Another ^He puzzle is the observed and appar- 
ently real dispersion in ^He abundances. Balser 
et al. (1994) present this dispersion as a ^He 
abundance gradient, with '^He decreasing towards 



the galactic center. This trend runs counter to 
intuition, if ^He is to be formed by low mass stars 
which are more prevalent towards the galactic 
center. Moreover, for stars having lifetimes this 
long, it is hard to see how this gradient could 
happen at all, as such stars will have traveled a 
significant fraction of disk in their lifetimes and 
so should be well spread throughout the Galaxy. 

A possible explanation for "what is wrong 
with ^He" was offered by Olive et al. (|1994D . 



They point out that the correlation of ^He with 
galactocentric distance might be better under- 
stood correlation with the mass of the H II 
region. Then the '^He in the higher mass regions 
will contain more eject a from high mass stars. In 
this scenario, one expects ^He to be lower wher- 
ever there is more ejecta from massive stars, i.e. 
more stellar processing. Thus the ^He abundance 
should go down towards the center of the galaxy, 
where the density and so the star formation rates 
are higher. Note that this could also answer the 
puzzle of the existence of the gradient despite the 
long lifetime of low mass stars. The point is that 
there is not a gradient in '^He production (by low 
mass stars), but instead there is a gradient in 
destruction by short-lived, high mass stars. 

For our chemical evolution models, we will 
want to use our apparatus to address the im- 
portant issue of the effect of low mass processing 
on ^He evolution, allowing for different possibili- 
ties as done by Vangioni-Flam, Olive, & Prantzos 
(1994). We will also include infall, as was done 
in Steigman & Tosi (|1992|) . 



2.5. Helium 4 

2.5.1. '^He Data 

Observations of ^He are well reviewed else- 



where, e.g. Pagel (|1993|) , Pagel et al. (|1992D , 
Walker et al. (p^9lD , Davidson & Kinman ( [1985D , 
Skillman & Kennicutt ( |1993 ); Skillman et al. 
(p9l ); and Copi et al. ([19951) . Briefly, accu- 
rate determinations of ^He are hard to come by. 
The best observations are of H II regions, where 
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He should be ionized and accurate abundances 
might result. Although there have been obser- 
vations of ^He in Galactic H II regions, systems 
with lower metallicity are expected to contain 
less contamination from stellar ejecta. The sites 
of choice have proven to be metal poor extra- 
galactic H II regions.^ 

The solar system abundance of ^He is Yq = 
0.274 ± 0.016 (Anders & Grevesse [19891) . ISM 
abundance determinations are few, and further- 
more, the ionization structure of gas regions 
causes ^He to be systematically underestimated. 



(Wilson & Rood 1994). Thus, we have no infor- 
mation on ^He in the ISM that is accurate enough 
to be an important constraint on chemical evo- 
lution models. 

2.5.2. ^He Chemical Evolution 

As the second most abundant nuclide in the 
universe, next to hydrogen, ^He comprises about 
a quarter of the universe's baryonic mass. It is 
also the only nuclide made in significant amounts 
in stars of a broad range of masses, with the dom- 
inant production in stars of middle masses (Ren- 
zini & Voh pJSl] ; Iben & Truran [l978| ) . Despite 
this nearly ubiquitous stellar synthesis, however, 
stars only make a small (< 10%) contribution to 
the total helium abundance (Hoyle and Tayler 



1964): most He is primordial. Nevertheless, to 



fully test BBN and moreover to use its power 
to constrain particle physics and/or test condi- 
tions in the early universe, one needs to know 
the abundance of ^He exceedingly well. 

Even in low metallicity H II regions, stellar 
pollution exists; the question is how best to de- 
termine it and so deduce the primordial He abun- 
dance. To this end Peimbert & Torres-Peimbert 
(|1974| , |1976| ) noted that He production should be 
correlated with heavier element production. In 



^Jakobsen et al. (L994) have recently reported detection of 
*He in a quasar absorption line system. The uncertainties 
in the measurement render it a qualitative confirmation 
of BBN, but with more accuracy this type of observation 
could be extremely powerful. 



H II regions one observes both He and CNO. Fol- 
lowing Peimbert & Torres-Peimbert, one deduces 
the primordial helium mass fraction from ex- 
trapolating the low metallicity end of the Y vs. 
Z plot, exploiting the relation 



Y 



Y +'^^ Z 



(17) 



Equation ( [ITD should be valid for sufficiently 
small Z (and dY/dZ a constant function of Z), 
for which this procedure should succeed in ex- 
trapolating Yp. The challenge for chemical evo- 
lution is to determine which Z are sufficiently 
small, and which "metals" C, N, and/or O are 
the best surrogates for Z. The challenge for ob- 
servers is to make precise enough measurements 
to allow for a meaningful extrapolation.^ 

Analysis using this technique has for the most 
part avoided the explicit use of detailed chemical 
evolution calculations. Instead, the approach has 
been an empirical one: Y and Z are measured 
for different regions; the Y ys. Z relation is plot- 
ted, fit, and extrapolated to get Yp = Y{Z = 0). 
Two complications to this procedure immedi- 
ately arise. First, one does not directly measure 
the metal fraction Z, and so one uses as a surro- 
gate either oxygen or nitrogen (and carbon — see 

|l98g| ). 



Sec- 



Steigman, Gallagher, &: Schramm 
ond, and moreover, one must adopt an expected 
Y vs. O (and N or C) relation to use in fitting 
the data and making the extrapolation. 

While everyone recognizes that the fits need 
not be linear, this is the simplest two parameter 
fit, and so this has been the first to be tried. 
A glance at the data (see figure suggests that 
this might not be a bad first guess. Indeed, the Y 
vs. 0/H relation is well described by a linear fit. 
In contrast, the case of nitrogen as a metallicity 



^The difficulties of measuring He and metal abundances 
to the desired accuracy (Fp good to the third decimal 
place!) cannot be overstated. For attempts to correct 
for some of th ese e rrors by culling points from t he data 
set see Pagel ( |l993| ), Pagel & Kazla uskas ( |l992[ ), Pagel, 
Simonson, Terlevich, & Edmunds, (1992), and Olive & 



Steigman (1995) 
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tracer has generated vigorous debate as to the 
appropriate fitting function. At issue is whether 
to choose a linear fit, in which the slope HsY/ 
is constant, or a nonlinear one, in which N oc K^, 
and thus the slope is not constant and has the 
possibility of being quite steep for small N, i.e. 



at early times (Fuller, Boyd, & Kalen [1991 ; but 
see Olive, Steigman, & Walker |99|). 

Put differently, the question regarding the 
y — N relation is, what is the N — O relation? 
The nonlinear fit for y vs. N arises by assuming 
that N production is proportional to the C and O 
abundance, leading to the "secondary" relation- 
ship N oc . We wish to know the relative size of 
the primary and secondary stellar contributions 
to N, an issue not yet resolved. Consequently, 
we will want to test effects of both primary and 
secondary N. 

Most ^He observations are determined from 
extragalactic H II regions. To completely under- 
stand these abundances would require separate 
evolutionary models for each observed galaxy. 
However, we will be fitting our model to the solar 
system, and so we should not necessarily expect 
a close fit to the extragalactic data. However, 
the few models that do exist, e.g. those of Math- 
ews Boyd & Fulle r (|l99^ ), and Balbes, Boyd, 
& Mathews ( |1993| ), are qualitatively similar to 
those studied here. 

Another outstanding problem in ^He evolu- 
tion is the theoretical reproduction of the slope 
Ay/AZ. The theoretical dY/dZ depends on how 
many high mass stars are included in the cal- 
culation (e.g. Maeder |1983| , |1992| , |1993| ; Brown 
& Bethe pl9^ ; and Prantzos pJ9^ ). Stellar 
yields for massive stars give a much lower ratio 
dY/dZ ~ 1.3 than that for intermediate and low 
mass stars, dY/dZ ~ 6. Averaged over a typical 
IMF with mass limits of 0.4 and 100 AIq gives 



dY/dZ - 2 



(18) 



(Pagel |T99|). 

One attempt to address this issue notes that 
many models for high mass stars (e.g. Weaver & 



Woosley 1993) do not include mass loss. For stars 
above ~ 30 Mq this can have a significant impact 
on all abundances but particularly ^He, as sug- 



gested by the results of Maeder (1992). Indeed, 
stars above ~ 30 Mq are thought to rapidly lose 
mass until reaching 3OM0; thus one need only 
compute stellar yields up to this mass. For our 
purposes, we will explicitly include the possibil- 
ity of mass loss in high mass stars. 

Finally, an unfortunate aspect of stellar evo- 
lution is that there is a great deal of uncertainty 
as to the gross behavior-let alone the nucleosyn- 
thesis yields-of stars in the 8-10 Mq mass range 
(further discussed below, § 3.1.1.| ). It has been 



at most (Maeder 1992, 1993). This theoretical 
slope is to be compared with an observed slope 
(calculated assuming Z oc O) of dY/dZ ~ 4 — 6 



suggested (Woosley & Weaver 1986 ) that these 
stars might dominantly produce helium. This 
could have a large effect on the net '^He produc- 
tion and so it is worth trying 8 — IOMq yields 
with very large yields helium, to see the sensitiv- 
ity to this range. 

3. The Suite of Models 

Our goal is to examine a broad selection of 
physically plausible scenarios, including at least 
the variations that have been commonly used in 
the literature. This approach is inspired by that 
of Tosi ( 1988| ), although her study made infall a 
central feature whereas we will emphasize other 
features as well. 

We classify model features as "global" in scope 
for particular to the light elements. We will des- 
ignate "global" model features to be: (1) those 
that affect all elements (e.g. choices of SFR or 
IMF) and the gas consumption; or (2) those that 
affect the yields of C, N, O, and Fe (CNOFe), 
i.e. the nuclides we will be following which are 
not light elements but serve as chronometers or 
tracers for the light elements. Upon identifying 
the set of global model variations, we combine 
all of these features independently to create a 
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large set of potential models for the solar neigh- 
borhood. By systematically investigating every 
possible combination of our chosen global fea- 
tures, we will be assured to properly estimate 
the full model variation inherent within the span 
of model options we have allowed. 

We demand that our results conform to the 
features of the solar neighborhood. By select- 
ing those models which reproduce the solar char- 
acteristics, we will develop a suite of acceptable 
global models. To each of these global models is 
then added a wide variety of light element evo- 
lutionary prescriptions. We then calculate the 
chemical evolution for every light element vari- 
ant, obtaining an inclusive measure of the model 
dependence of light element evolution. 

3.1. Selection of Global Model Features 



In her study, Tosi (|l98j ) found that all aspects 
of chemical evolution models affect the absolute 
abundances. However, she found that the choice 
of IMF, and of course the selection of nucleosyn- 
thesis yields, had the most profound influence, 
affecting the abundance ratios as well as the ab- 
solute abundances. Therefore we will want to 
be particularly thorough in our consideration of 
these model features. Our global model proper- 
ties are summarized in table ||. 

3.1.1. Nucleosynthesis Yields 

The tabulated nucleosynthesis yields are those 
of low to intermediate mass stars, with masses 
~ O.8M0 to ~ 8M0 , and high mass stars with 
masses M > \2Mq . The former are understood 
to be stars whose lives progress to the point of 
planetary nebula ejection and becoming carbon- 
oxygen white dwarfs; the latter, high mass stars 
end their lives explosively as type II supernovae. 
Standard results for this study are those of Ren- 



zini & Voli ( 1981 ) for intermediate mass stars, 
and Weaver & Woosley ( 1993| ) for high mass 
stars. 

The uncertainties in the yields pointed out in 



these references. The Renzini &: Voli (1981) re- 
sults, for example, were found to depend on the 
adopted parameterization of convection in the 
form of the mixing length. However, this level 
of detail is not a major source of uncertainty for 
our purposes. It would be desirable to compare 
results for different stellar models and yields; un- 
fortunately this is not possible for the low mass 



stars, as Renzini & Voli (1981) are the only cur- 
rent models with detailed results reported. In 
addition, for the massive stars, several groups 
have published massive star yields; however these 
models have very different degrees of emphasis 
on their nuclear yields, with Woosley & Weaver 
(1993) presenting particularly meticulous results. 
As motivated in §12.5. 1, we allow for a wide 



range of N evolution, including both primary and 
secondary sources for N. We thus will follow the 
approach of Mathews, Boyd, & Fuller (|1993 ) in 
writing the mass fraction of ejected N, X'^J, as 



a + [3 



(19) 



where Xq is the usual N yield, Xc{t) is the calcu- 
lated C mass fraction and X^ is the solar carbon 



mass fraction from Anders & Grevesse (1989). 
The constants a and /? are chosen as follows: 

(a,/3) = (1,0) ; (0.5,1) ; (0,1) ; (0,2) (20) 

where the first option is the standard case of 
purely primary N, the second is a mixed case 
of some primary and some secondary N, and the 
final two cases are variation on purely secondary 
N. Note that our procedure differs somewhat 



from that of Mathews, Boyd, & Fuller (|1993D 
who used O as the seed nucleus for secondary 
N, instead of C. Carbon is the more appropriate 



choice (Audouze, Lequeux, & Vigroux 1975 ; Vi- 
groux, Audouze, & Lequeux |1976 ; and Dearborn, 
Tinsley, & Schramm 1978| ), and we have used it. 



3.1.2. Initial Mass Function 

The selection of an initial mass function is a 
crucial feature of any chemical evolution model. 
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Unfortunately, there is no convincing physical 
theory of star formation, and so the nature of the 
IMF and the star formation rate over the history 
of the Galaxy are poorly understood. It is un- 
clear, for example, whether the IMF has changed 
with time, what its upper and lower mass lim- 
its are (and have been), and whether it depends 
upon the composition of the gas which will be- 
come stars. Furthermore, it is very difficult to 
untangle the behavior of the IMF and the SFR, 
which could in fact be inseparable and better left 
in the form of a creation function, as in eq. (^). 

Lacking any good theoretical or observational 
guidance as to these questions, we will adopt 
eq. @, i.e. we will assume that the IMF is con- 
stant in time. We will, however, allow for it 
to take different forms. We will try the classic 
Salpeter ( [1955| ) function, 



m oc m 



(21) 



and we will allow the index ("slope") x to vary, 
trying x = 1,1.35, and 1.7. We will also try 
an IMF (p derived from the observed, present- 
day mass function (PDMF, 4>^'^), as given by, 
e.g. Scalo ( |1986| ). The relationship between the 
two depends on whether the stellar lifetime is 
short enough that some of the stars of mass m 
have died; specifically. 



[m 



mj 



ct>^H^)/ jt'dt bit), 



> to 

m ^ to 



(22) 

where the age of the Galaxy is to and h{t) = 
'i/'(t)/(V')- We remind the reader that this proce- 
dure is only practicable when the star formation 
rate "0 is a given, explicit function of time,^ be- 
cause of the role of ip in eq- (p2|). 

We also will try different mass limits to the 
IMF. Here we are guided by common choices that 



In fact, other SFRs can accommodate the Scalo IMF 
through an iterative process (Mathews, Bazan, & Cowan 
1992); such a procedure would however be very time con- 
suming to use with our many models, and is not consid- 
ered here. 



have appeared in the literature. We will wish 
to have a broad, medium, and narrow range in 
mass. We choose limis of 

{mi,mu) = (0.2, 100); (0.1, 60); (0.4, 30) (23) 

which we feel to be representative. 

3.1.3. Star Formation Rate 

As the physics of star formation is ill under- 
stood, the form of the SFR is not well set by 
theory, although there have been some attempts 
to do so. As with the IMF, there are many vari- 
ants, but we will select what we feel to be typical 
choices: 



ftot 



l^CTtot 



V^'toty 



a exp(— t/r) 



(24) 

(25) 

(26) 
(27) 



where we pick the r = 7.5, 15 Gyr as the 
timescale for the exponential case. Note that 
these always decrease for a closed box (ctot = 
const), but that for infall models with accreting 
gas, these will show a rise at early times. 

3.1.4. Infall 

There is very little known about the existence, 
nature, and composition of infall and/or outflow 
from our Galaxy. Evidence for infall is perhaps 
given by the observation of the High and Very 
High Velocity clouds. However, there is certainly 
a need for infalling matter at early times, if one 
views infall to disk as just the outflow from the 
halo as it collapsed. In this case, one imagines 
a metal-poor (i.e. BBN composition) infall on a 
short (~ 2 Gyr) timescale. Another justification 
for infall is more pragmatic: it provides a pos- 
sible solution to the G-dwarf problem (discussed 
in §111). 

We will adopt models with and without infall, 
which we parameterize according to the widely 
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used form f{t) = /o exp(— t/rinf). We will vary 
the infall strength /o, and the timescale Tinf. 
Namely, we will try Tjnf = 2, 4 Gyr, and we will 
use different /o such that infall contributes 50% 
and 99.9% of the total disk mass today. We will 
also consider "closed box" models with no infall, 
i.e. for which /o = 0. 

3.1.5. Supernova Rates 

We will assume that all supernovae arise ei- 
ther from the core collapse of massive stars (type 
II) or from accretion onto a CO white dwarf 
in a binary system (type la). All stars above 
~ 12Mq are thought to become type II super- 
novae, whereas the CO progenitors of type la su- 
pernovae are arise from binary systems of inter- 
mediate (3 — 8 Mq) mass. The rates for these 



putting 



E., 



la 



(29) 



events are given by (Matteucci & Greggio 1986 ) 
and are of the form 



(28) 



where J- \s a. functional of the IMF and SFR, as 
well as /, the primary-to-secondary mass distri- 
bution in binaries (and perhaps also the initial 
separation; e.g., Mathews et al. 1993| ). 

The amplitude of the type la rate is controlled 
by the dimensionless parameter A which is the 
probability for binary systems in the appropriate 
mass range to undergo supernova events. This 
number is not well understood theoretically and 
in practice is adjusted to reproduce the present 
ratio of type II to type la events ~ 1 ; this ratio is 
proportional to A but is much larger due to the 
larger number of low-mass, potential type la pro- 
genitors. We will examine the sensitivity of our 
results to this parameter, by choosing two values 
for A: a high value of A = 0.05 near the value rec- 
ommended by Matteucci & Greggio ( 1986|) , and 
a low value A = 0.007 used by Timmes, Woosley, 
& Weaver (|l995| ). 

The type II nucleosynthesis yields are those 
discussed in § |3.1.1. , and are directly included in 
eq. (^) . Type la yields of nuclide i are added by 



with M^^ the mass ejected in i in a type la event. 
We use the yields for the W7 model of Thiele- 
mann, Nomoto, & Yokoi ( [L986| ). 

3.2. Constraints on the Global Models 

We first constrain the global models to fit the 
solar neighborhood properties within fairly gen- 
erous tolerances, excluding for the moment any 
consideration of the acceptability of the global 
models' light element yields. The available con- 
straints on the gross features of the models are 
(Tg, /ii, and the CNOFe abundances. 

Note that we want to be not too restrictive 
on our models, as the data, even the solar sys- 
tem abundances, are likely to be subject to sys- 
tematic errors. For example, Olive & Schramm 



( 1981 ) suggest that the solar system may have 
been formed from a young OB association; as a 
result its elemental and isotopic content would 
be a biased indicator of the larger solar neigh- 
borhood. Indeed, it seems unwise to dwell on 
one constraint at the exclusion of the rest. Con- 
sequently, we wish to be generous in choosing 
models with which to check chemical evolution 
uncertainties in the light elements. 

We will demand that all models fit the lo- 
cal gas mass fraction, which is determined to be 
^ = 0.1 — 0.2 (Rana 1991| ). Specifically we re- 
quire /i = 0.13 — 0.17. This fixes the overall mass 
consumption. Given that the a central feature 
explicitly built into chemical evolution models is 
the mass consumption and the conservation of 
overall mass, this is a most important constraint. 
We therefore have tuned all of the SFR normal- 
izations to satisfy this constraint. 

As for the quantitative means of evaluating a 
model, the analysis of Tosi ( 1988| ) is interest- 
ing but perhaps premature given that the data 
is likely to be fraught with systematic errors. In- 
stead, we will make the generous (and implemen- 
tationally simple) demand that the solar CNOFe 



14 



abundances within a ±0.4 dex^ range. This is to 
be compared with the factor of 2 that Timmes, 
Woosley, &: Weaver ( |1995| ) allowed themselves in 
a model tuned to give good results. 

We will also constrain some of our models with 
observed G-dwarf distribution. G-dwarf stars are 
long lived, with ages comparable to the Galactic 
age now. Consequently, their provide a record 
of the integrated star formation history of the 
galaxy; the metallicity distribution of these stars 
is as seen in Figure ^. It has long been known 
that the simplest models of galactic chemical evo- 
lution are unable to reproduce the G-dwarf distri- 
bution. In particular, such models cannot avoid 
overproduction of low metallicity stars, whereas 
the data shows a dearth of G-dwarfs below a 
metallicity of [Fe/H]. 

One of the solutions to this classic "G-dwarf 
problem" is the introduction of infall. In models 
with infall, low metallicity stars are still formed, 
but as material is added subsequently, the frac- 
tion of low metallicity G-dwarfs becomes accept- 
ably small. We will thus apply a G-dwarf con- 
straint to our models, but only on models with 
infall. For these we will demand a solution to 
the most emphasized part of the G-dwarf prob- 
lem, namely the need for a small number of low 
metallicity G-dwarfs. We will ask that our mod- 
els not exceed the observed fraction of such stars 
at [Fe/H] = -0.8 by more than a factor of 0.4 
dex ~ 2.5; this will prove to be a restrictive con- 
dition. We will not make the stronger demand 
that the models fit the full G-dwarf distribution 
to a high accuracy, as this is not done even in 
models which allow themselves tuning which we 
lack. 

We will not impose the G-dwarf constraint on 
models without infall, as it has long been thought 



(e.g. Tinsley |1980| ; Pagel |1989|) that such closed 
box models are unable to meet this constraint 



^dex = base 10 logarithmic units; i.e. a logarithmic vari- 
ation by ±0.4 dex is a linear variation by a factor of 



10" 



2.5 



without further modification. On the other hand, 
Mathews &: Bazan ( 1990| ) show that the G-dwarf 
problem is somewhat mollified by the inclusion of 
metallicity dependent stellar ages, and Sommer- 
Larsen ( 199 1] ) has noted a selection bias in fa- 
vor of high metallicity dwarfs. At any rate, we 
keep the closed box models since their ubiquitous 
presence in the literature demands that one ex- 
amine their results, and again, we want to err on 
the side of inclusiveness in our selection of viable 
models. 

3.3. Light Element Model Features 

Having assembled a set of possible models and 
determined constraints on their gross character- 
istics, we now turn to the model features we will 
use to encompass the possibilities for light ele- 
ment evolution. For the most part, these intro- 
duce small enough changes in the CNOFe abun- 
dances that each model's CNOFe prediction will 
remain valid despite the modification of light ele- 
ment evolutionary features. An exception to this 
rule comes from mass loss. 

The evolution of all light elements is affected 
by their initial BBN abundance as set by the 
choice of r/. We use three different values of 
r], guided by the analysis of Copi, Schramm, & 
Turner 1995 (cf. eq. |l]). We choose central, high, 
and low values r/ = (2,4,6) x lO"^'^; these values 
bracket the Copi, Schramm, & Turner's "sensi- 
ble" limit. The low limit is intentionally cho- 
sen to be slightly below their recommended value 
(but above their "extreme" lower limit), in order 
to allow comment on recent speculation of the 
possibility of a low r]. 

In addition to systematically varying 77, we will 
include the following features for each element: 

3.3.1. D 

As we have discussed, D evolution is uniquely 
simple because there are no stellar sources of D; 
the range in D results thus stems only from the 
choice of rj and from global model features. 
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3.3.2. ^He 

For high mass '^He yields we used the results 



of Weaver &; Woosley ( 1994 ), which span a grid 
of masses from ~ 1 1 — 40 Mq and very complete 
range of metallicity from Z = to Z = Zq . The 
metallicity dependence of these models is strong 
and thus a crucial feature to include; these re- 
sults update the work of Dearborn, Schramm, &: 
Steigman ( |1986| ). 

As noted in § ^.4j , '^He evolution hinges on the 
effect of low mass stars on ^He. To investigate 
this effect we have chosen an approach similar to 



that of Vangioni-Flam, Olive, Sz Prantzos (1994), 



as well as Olive et al. ( 1994 ). Namely, we specify 
the ^He survival fraction at m = 1, 2, and 
3 Mq. We choose to allow for: (1) modest ^He 
destruction (Vangioni-Flam, Olive, & Prantzos 
found they must); (2) "break even," with all '^He 
surviving but no net production; and (3) modest 
^He production. Specifically, the three ranges we 
choose are 



53(1M0,2M0,3M0 



(1,0.7,0.7) 
(1,1,1) (30) 
(1.33,1.33,1.33) 



We chose the production to be small (i.e. some- 
what less that the Iben &: Truran ( |197^ ) yields, 
corrected to include initial D, would suggest) in 
order to be conservative. It is clear, as pointed 
out in Olive et al. ( |1994| ) that the yields at face 
value do not work, assuming the conventional 
^He high mass yields. In an effort to allow for 
uncertainties in the Iben & Truran ( |1978| ) calcu- 
lation, we wish to see if even a modest production 
can work. 

3.3.3. ^He 



As discussed in § 2^ , calculations of ^He pro- 
duction must allow for its production by stars 
having a wide range of masses; because of this 
broad range, ^He is particularly sensitive to sys- 
tematic problems in nucleosynthesis yields over 
the whole mass range. One such problems is that 



tabulations of stellar nucleosynthesis yields usu- 
ally do not include results for stars in the range 
of roughly 8 — 12 Mq. There is not a firm under- 
standing of the basic nature of these stars' final 
fate; whether they become white dwarfs or ex- 
plode as supernovae is a delicate and unresolved 
question. The nucleosynthesis results for these 
stars are therefore quite uncertain, and indeed 
possibly quite different from those of stars ei- 
ther more or less massive. This is a regrettable 
state of affairs, as a typical IMF places a good 
deal of weight on the yields from stars in this 
range. Thus whatever assumptions one adopts 
about these yields can prove important for chem- 
ical evolution, and so we will examine two possi- 
ble scenarios. 

Traditionally, one obtains yields in this mass 
range by interpolating the results from the closest 
tabulated masses. Such an interpolation can be 
viewed as hedging one's bets about the accuracy 
of estimates for the onset of core collapse, and we 
will include this case; however, it has also been 
suggested that stars in this range produce ^He al- 



most exclusively (Woosley & Weaver 1986 ). We 
therefore mock up this behavior by assuming that 
all of the ejecta from these stars is in this form. 

Other important model uncertainties regard- 
ing '^He production arise due to the omission 
of mass loss in high mass stellar models Con- 
sequently, we use the stated Woosley & Weaver 



(|1993|) high mass yields, which do not include the 
effects of mass loss, but we also characterize the 



effects of mass loss by adding Maeder's (|1992D 
calculation of these. Note that while Maeder also 
calculated supernova yields which included mass 
loss, the models were far less detailed and thor- 
ough than those of Woosley & Weaver; thus we 
choose to import only the mass loss addition. 

4. Results 
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4.1. Effect of the Chemical Evolution Model 
Features 

The effect of our global model features are for 
the most part well documented in the literature; 
we will discuss only the highlights. For a more 
complete discussion we refer the reader to, e.g., 
Talbot & Arnett (|197lD , Tinsley (|198q ), and Tosi 
( |1988| ). Note that unless otherwise indicated, all 
of the figures in this section are for models which 
vary only the parameters discussed while keeping 
the others constant. 

A sample of SFR behaviors appears in Fig- 
ure ^; Figure § (a) displays the case of a closed 
box model. We see that for SFRs of the form 
tp = z^fTtot/^" oc cjg", there are progressively 
higher initial bursts and more rapid declines for 
progressively higher n. The presence of infall 
changes the behavior of the SFRs whose form 
depends upon the total and/or gas masses. As 
seen in Figure |2| (b), for ip = uatot the SFR is 
monotonically increasing with time as the disk 
mass grows through infall, in contrast to the de- 
caying SFR forms in a closed box. For ■0 = vag 
(Figure ^ (c)), in the case of infall the SFR ini- 
tially rises to a peak, but then falls off as the gas 
is converted to stars and remnants. 

The total and gas mass evolution is illustrated 
in Figure ^. Our SFR boundary conditions and 
infall prescription fix the fitot plot as well as the 
endpoints of the cTg plot. Note that the spread in 
(Tg at ti = 15 Gyr is an indication of this spread 
in Hi, which arises due to slight differences in the 
SFR tuning for different models. The effect of in- 
fall on the gas and total mass is clear from Figure 
^. Without infall, the total mass is of course con- 
stant, and the gas mass monotonically decreases. 
With infall, the total mass rises rapidly, with the 
gas mass initially following, then reaching a peak 
and turning over. 

The relation between iron abundance and time 
(the "age-met allicity relation") is an important 
one: iron is readily observable whereas t is not, 
and the maximal binding energy of Fe ensures 



that its abundance will increase monotonically 
with time. Thus the iron abundance is tradition- 
ally used as a chronometer (although its origin in 
both type la and type 11 supernovae makes it a 
somewhat problematic one, with O perhaps be- 
ing a more suitable choice. See Wheeler, Sneden, 
& Truran |1989| ). 

In Figure |^ we plot [Fe/H] versus t for all of 
our SFR options, with data from Edvardsson et 
al. ( 1993| ). In general, we see that the Fe rises 
very quickly, with most of the production hap- 
pening at very early times and very little new Fe 
being added today. In comparing the curves in 
Figure |^ to the data it is immediately clear that 
all of the models are viable. There is also a great 
deal of scatter in the data, indeed, as much as the 



overall trend. Edvarsson et al. (1993) emphasize 
that this scatter is real, and is not well explained 
in most chemical evolution models which invoke, 
as we have, the assumption of instantaneous mix- 
ing of stellar ejecta at a given epoch. We can only 
expect to reproduce the average trends in the ob- 
servations, and one should recall that individual 
stars can show considerable excursions from the 
average trend. 

The evolution of deuterium provides a good 
example of the effect of different IMFs; the D 
history for three different IMFs is displayed in 
Figure ^ (a). The three choices have different 
return fractions, with R = 0.32, 0.42, and 0.57. 
As expected (cf. eq. |^), the D evolution is very 
sensitive to this parameter, with the overall de- 
pletion at the solar system (i.e. at = 10.4 Gyr) 
varying from a factor of about 1.4 to a factor of 
2.3. As we will see, this spread accounts for a 
large portion of the model uncertainty in the D 
evolution (for a fixed rj). Note that the largest 



D depletion comes from using the Scalo (1986) 
IMF, which is derived from the present-day mass 
function. This IMF is bimodal and so has a lot 
of power in the heaviest stars which eject most 
of their mass; consequently, this yields a large 
return fraction and considerable D depletion. 

In Figure |5| (b) we show ^He and D-|-^He for 
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the same models. The tradeoff between D and 
^He is clearly seen, with ^He being progressively 
higher in models having progressively lower D. 
Indeed, despite the large variation in the individ- 
ual D and '^He cases for the three IMF choices, 
we see that all show a very similar — and nearly 
constant — evolution of D+^He. In Figure ^ (c) 
we show the ^He evolution for the same models. 
These models show a correlation between the D 
and ^He evolution, as the models with larger D 
depletion also have larger ^He production. We 



will investigate this relationship in detail in §4.3. 

Deuterium together with ^He provides a good 
barometer of the effects of infall on the evolution 
of an isotope. The infall of primordial material 
enriches the depleted Galactic D, while it dilutes 
the increasing Galactic ^He. This behavior is ev- 
ident in Figure |6| which shows the slower decline 
of D in the infall models, mirrored by the inhib- 
ited growth of ^He. Note, however, that in the 
case of D the spread (at the time of the formation 
of the solar system) due to different treatments 
of infall is only about 15%, much less than that 
caused by different IMF choices. 

The different behaviors of nitrogen in our 
models are illustrated in Figure ^ (a). We focus 
on N behavior at the same low metallicities which 
are observed extragalactic H II regions, and we 



include data from Pagel et al. ( 1992 ) and Skill- 
man et al. (1994). While there is a large scatter 
in the data, the trend seems better fit by the 
models with mixed and the purely primary N. 
Note in particular that the mixed model does 
a surprisingly good job at the low metallicities. 
While the fit is to somewhat an accident of our 
particular choices for the primary /secondary N 
mix, it is amusing that the lowest metallicity 
region — which dominates the l^fit — comes out as 
well as it does. 

The effect of secondary N on the N relation 
is shown in Figure ^ (b), with the Y—0 relation 
in Figure (c). For these plots we have chosen 
models with elevated ^He production, as this pro- 
vides a better fit to the data, ( ^4.4.| ). The models 



with pure secondary N do indeed have a sharp in- 
crease at low metallicities; this arises from helium 
production in concert with suppressed N produc- 
tion due to the dearth of its seed nucleus C. Note, 
however, that the low metallicity y— N relation 
is much more gentle for the case of partial pri- 
mary and partial secondary N, the case that best 
fits the observations. The small departure from 
linearity in this case suggests that linear Y— N 
fits might not be in much error (§ |4.4.| ). 

4.2. Global Constraints and Model Un- 
certainties 

Having sketched some of the effects of the 
global model features, we now explore the range 
of abundances calculated for the initial suite of 
candidate global models and for the subsequent 
set of light element models. In finding the range 
of light element predictions we obtain a rough 
quantitative estimate for the model uncertain- 
ties inherent in the chemical evolution arguments 
used to deduce the primordial abundances of 
these elements. 

We created the suite of candidate global mod- 
els by making every allowed combination of fea- 
tures shown in table this produces 1184 mod- 
els. We ran all of these models, using one partic- 
ular choice of light element model features. Re- 
sults of this run are summarized in Figure |^ (a) , 
in which we plot the abundances for each ele- 
ment, calculated at the birth of the solar system 
{tQ, 4.6 Gyr ago). Abundances are expressed in 
terms of the logarithm with respect to the ob- 
served solar abundance A/H5^ 



obs. 





[A/H] = log 



A/H(te) 



(31) 



(solar data is from Geiss ( |1993|) fo r D and ^He, 
and from Anders & Grevesse ( |1989| ) for all other 
elements). Note the wide range of all calculated 
abundances, including those for the light ele- 
ments. Indeed, the only element without large 
excursions is ^He, since it is the only element 
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whose solar system abundance is mostly primor- 
dial and so on this scale does not show large ef- 
fects of uncertainty due to chemical evolution. 

We now impose the constraint that our models 
reproduce the solar abundances of CNOFe within 
±0.4 dex, i.e. within a factor ^ 2.5. While this 
may seem generous, recall that we have not tuned 
our models to fit any elements in particular, yet 
even models that do so (e.g. Timmes, Woosley, 

Weaver 1995| ) allow themselves a factor of 2. 
Also, the scatter in the age-metallicity relation 
suggests that there is a large spread of metallic- 
ities at a given time. We must be aware of the 
possibility that the solar data will not necessar- 
ily reflect the average at the solar neighborhood, 
so we do not want to exclude potentially good 
models on the basis of a single element. 

When we impose this constraint on the can- 
didate global models, we reduce the field to 267 
models. A plot of the solar abundances for these 
models appears in Figure ^ (b). The light el- 
ement uncertainties are seen to be greatly re- 
duced, while still being considerable. The spread 
in D and ^He is of particular interest, as we have 
not added any model variance in the D or ^He 
evolution (e.g. different rj or different low mass 
^He stellar processing). Thus the range abun- 
dances of D and ^He reflect only the variance due 
to the global models. Quantitatively, we have a 
variance in D of a factor of 2.0. As we have noted, 
this is very close to the range due solely to the 
IMF, which we see to be the largest source of D 
variation. We also find ^He to vary by a factor 
of 1.9, and AY to vary from 0.018 to 0.046. 

As discussed in 93.1.4J, we also constrain the 



models with infall to fit the G-dwarf distribu- 
tion, again with a loose tolerance. We focus on 
the ability the models to avoid a large number of 
low metallicity G-dwarfs, comparing our results 



to the data of Rana ( 1991 ) . A plot of the distri- 
butions for one model that passes and one that 
fails can be found in Figure P. As discussed in 



the number of models to 114, with 28 of these 
being infall models.^ 

We now have a set of 114 chemical models 
which fit the solar system behavior within the 
tolerances we have set. To investigate the range 
of light element evolution allowed by these mod- 
els, we run each of them for 18 different varieties 
of light element evolution models. Note that this 
is far less than the possible number of allowed 
combinations of options for light element evolu- 
tion. However, to run such a number of models 
would not only be tremendously time consuming, 
but it would also be redundant. We do not ex- 
pect, for example, the various LiBeB yields in the 
cosmic ray models to depend on rj in an impor- 
tant way compared to the other large uncertain- 
ties. Thus we have reduced the number of light 
element variations to only those 18 most likely 
to give interesting results. In total, therefore, we 
have 18 light element models x 114 global models 
= 2052 overall models to run. 

Results for the 2052 models appear in Figure 
10. Notice that a problem has developed in the 
CNO abundances — some light element evolution 
features have affected them. The problem arises 
due to the extra CNO contribution from massive 
stars with mass loss. For consistency with our 
protocol of constraints, we again constrain the 
models to give CNOFe within ±0.4 dex of solar. 

This reconstraint reduces the number of mod- 
els by about 25% to 1460. Results appear in 
Figure 11 (a). Note that even for this smaller 
number of models, the range in the light element 
abundances remains very similar, although with 
fewer models the extreme variations are less pop- 



ulated. The models depicted in Figure 11 (a) will 



3.1.4., we do not apply this constraint to models 
without infall. With the constraint, we reduce 



^ The only single features clearly important in determin- 
ing whether a given model would pass or fail were the 
IMF slope and mass limits. We found no models with 
X = 2 survived; for x = 2.35, no models with the widest 
mass limits passed; for x = 2.7, mostly the models with 
the smallest mass limits passed. Aside from the IMF, no 
other model feature was prominent in either its presence 
or absence from the passing models. 
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be the basis for the analysis in the sections be- 
low, except where otherwise noted; we will refer 
to these as the "final set" of models. 

The source of the gap in the D abundances 
in Figure ^ (a) becomes clear when we plot the 
results for a particular tj. As Figure |1^ (b) shows, 
the D gap arises due to the very high initial D 
at r/ = 2 X lO"^*'. It is also clear that it is hard 
for our models to reduce such a high initial D to 
its presolar abundance, a point we will explore in 
more detail in the next section. Note that ^He 
also shows sensitivity to its initial abundance. 

The spread in the light element abundances, 
particularly the spread beyond that in the "con- 
trol" group CNOFe, gives a sense of the spread 
due to model dependence. Note, however, that 
one must be careful in interpreting these limits. 
While we have tried a broad sample of models, 
the lack of a firm theory of star formation leaves 
open the possibility that one might construct 
other ones. Moreover, recall that we have not 
used chemical evolution to get limits on the pri- 
mordial abundances. Instead, we have assumed 
BBN yields and seen how well these fit the solar 
system. What one really would like to know is 
not this but rather the opposite case. Namely, 
one would construct models with variable pri- 
mordial light element abundances (i.e. not nec- 
essarily given by BBN values at a given 77), and 
constrain these models to reproduce the solar 
system values for the light elements and every- 
thing else. The resulting spread in the primor- 
dial abundances would be an indication of the 
model uncertainties incurred when setting lim- 
its on r]. This procedure suffers from an exclu- 
sive reliance on the veracity of chemical evolu- 
tion models. Furthermore, we may estimate the 
model ranges such a calculation would give, as 
we now discuss. 

Table | summarizes the model variations in 
the light element production for each 77. We see 
that, for all rj, D varies almost exactly by a factor 
of 2 and Ay varies by 0.063 units. The vari- 
ation in D+^He depends somewhat on rj, but 



is small, varying by < 20% from ~ 2.8. The 
r/-independence of these quantities are evidence 
that the range model variations does not depend 
on the light element initial values. This gives 
us confidence that, were we to do a light el- 
ement evolution study by allowing the primor- 
dial values to float (as described in the preceding 
paragraph), we would find similar ranges. That 
is, we expect the reduced quantities in table ^, 
namely Dmax 

/B^i^, (D+3He)„,ax/(D+3He)^in, 
and Ytnax — ^mim to be generic estimates for 
the chemical evolution model ranges of these el- 
ements. Of course other selections of model fea- 
tures will give different numbers, but these re- 
sults at least give a quantitative estimate for the 
model dependence within a large scale of models, 
and serve as a benchmark for others. 

We emphasize that the D variation is due 
solely to variation of global model features, in 
particular the IMF and also infall. Thus to re- 
duce the variation in D requires better chemical 
evolution modeling clS Si whole. By contrast, the 
span in the '^He and ^He ranges are due in part 
to the effects of different stellar processing fea- 
tures for these elements. Improvements in these 
features could reduce the ranges given, indepen- 
dently of improvements in the global chemical 
evolution framework. 

It is encouraging to note that the variations 
in D, '^He, and ^He, while large, are significantly 
smaller than the 0.8 dex (i.e. factor of 6.3) over- 
all variations allowed for CNOFe in our global 
model constraints. This implies that the light el- 
ements through ^He have less model uncertainty 
than do the heavy elements. However, the re- 
lationship between the heavy and light element 
uncertainties is not necessarily a simple one. 

To get a feel for the connection between the 
heavy and light element model variations, we ex- 
amine the effect of tightening the fiducial toler- 
ances in the heavy elements. Specifically, we now 
require that our models reproduce solar CNOFe 
to within ±0.3 dex, i.e. to within a factor of 



2. Results appear in Figure 12. The reduction 
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in the number of models is drastic, going from 
1460 models at ±0.4 dex tolerance to 479 mod- 
els at ±0.3 dex tolerance, a 66% reduction in 
models from only a 20% tightening of the con- 
straint. One can understand this large reduction 
by noting that our heavy element yields generi- 
cally are scattered around the solar levels. For 
example ^^C is produced on "average" above the 
solar level, while ^^O is produced in a large range 
centered below solar; in other words, the 0/C ra- 
tio is typically subsolar. Thus our models have 
some range of error in reproducing O/Cq, and 
will never have all abundances any more accurate 
than this range. Indeed, demanding simultane- 
ous fitting of both abundances for progressively 
smaller tolerances would eventually eliminate all 
models when the constraints become too strin- 
gent given the heavy element yields we have used. 

Despite the large reduction in number of mod- 
els which survive the stronger heavy element con- 
straint, the spreads in the light element abun- 
dances remain roughly the same. This reinforces 
the notion that the CNOFe spreads are to some 
degree decoupled with the light element spreads. 
More quantitatively, we find for these models 
that D varies by about a factor of 1.76 and ^He 
by a factor of 3.1, values closer to the range for 
models with the ±0.4 dex CNOFe tolerance than 
the 20% tightening would naively allow. In other 
words, by making the variation in CNOFe smaller 
by 20%, we do not reduce the variation in D and 
^He by as large a factor. Thus the model varia- 
tion in the heavy and the light elements are not 
related by a simple linear scaling, but instead the 
relationship between the two seems to be more 
subtile. 

4.3. Correlations in Light Element Evo- 
lution 

We now study in detail the results from the 
set of 1460 models comprising the "final set" hav- 
ing all allowed light element variation while also 
fitting CNOFe. In particular, we will be inter- 
ested looking at D, ^He, and ^He (pre)solar and 



ISM predictions simultaneously to see better the 
correlations in their evolution and to show the 
interplay of the observational constraints. 



In Figure |1^ we plot the presolar D abun- 
dance versus the presolar ^He abundance for each 
model. One notices immediately that the choice 
of r] (different shaped points) is the largest ef- 
fect controlling the presolar abundance of both 
elements. Furthermore, the correlation between 
D and ^He is evident: the negative slope in the 
lines traced out in the Figure points out the well- 
known tradeoff that D destruction leads to ^He 
production. The steepness of these lines (which 
radiate out from the primordial values at each 
r]), depends on the degree of low mass destruc- 
tion or production of '^He. The highest of the 
three "prongs" for each r] comes from the case of 
^He production in low mass stars. The middle 
prong is the "break even" case, and the lowest 
one from net ^He destruction. We clearly see 
that the ^He presolar yields are very sensitive 
to the low mass stellar behavior, which we have 
only varied by 30% either way from break even. 
Our models thus confirm the results of Vangioni- 
Flam, Olive, k Prantzos ( |1994| ) and Olive et al. 
(1994), and extend them to include a very wide 
range of models. 



Figure |13| also includes lines marking the 2 — a 
variation in the observed presolar D and ^He 
abundances, as calculated by Geiss ( [1993| ). The 
Figure provocatively shows that none of the low- 
r], i.e. high initial D and '^He, models are able 
to fit the presolar data to within this tolerance. 
Indeed, the majority of the moderate-r/ models 
also seem to fail. At face value, this might sug- 
gest that low values of tj are ruled out and that 
higher ones are favored. The question is whether 
this strong conclusion holds up under scrutiny. 

As we have discussed in the previous section, 
it is a subtle matter to use our models to con- 
strain the possible BBN abundances. We have 
not included all possible chemical evolution mod- 
els (e.g. IMF and SFR schemes) and one cannot 
conceivably do this, given the freedom one has 
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to adjust each. Also, recall the observed scat- 
ter in Fe abundances for stars which are presum- 
ably coeval. If the presolar D is a low varia- 
tion to the mean we have calculated, we would 
be incorrectly constraining our results. Finally, 
as we have noted, significant uncertainties and 
assumptions underlie the entire chemical evolu- 
tion framework we have adopted (and variants 
thereof). Even if one could show that if no such 
models work to explain the low-?7 D and ^He, 
such a conclusion would not be founded upon a 
first principles calculation, since none exists. 

Nevertheless, it is not terribly surprising that 
7/ = 2 X 10^^'^ fails for our study. Indeed, as we 
have noted, this value is outside of the "sensible" 
range in rj given by Copi, Schramm, &: Turner 
( |1995| ), although it is within the "extreme" range. 
Also, while chemical evolution models are uncer- 
tain, as we have argued, the D evolution is per- 
haps the most certain of any nuclide considered. 
Thus there is no obvious reason to disbelieve the 
D results, particularly when other abundances 
can be fit reasonably well. 

Furthermore, while chemical evolution is known 
to be uncertain, our results as presented help one 
to quantify the size of the model dependence. 
At 7? = 2 X 10"^'', the D abundance varies by 
a factor of ~ 2; but the lowest D abundance for 
this r] is still a factor of ~ 1.5 away from upper 
limit to the solar system observation (and with 
an undesirable ^He abundance at that). It is in- 
correct to think of the spread in our results as 
some sort of distribution around a mean model, 
as we have not just performed smooth parame- 
ter variations, but we have also changed whole 
prescriptions for, e.g., the SFR. Thus it is hard 
to say more than that the presolar value of D 
for BBN with r/ = 2 x 10^^'^ is not reachable in 
our models and seems hard to reach by models 
similar to ours. 

While presolar D is a strong constraint on 
models with high r], we see that presolar ^He is 
a very strong constraint for all 7]. Because ^He 
survives processing of both D and ^He, its abun- 



dances grows with time.^ However, the presolar 
^He is not much larger than its primordial level, 
and so demands that models produce very little 
^He. We see that models with high rj, i.e. with 
the lowest initial ^He, are the only ones that can 
fit the presolar constraint. Within these, the dif- 
ferent ^He processing becomes important, with 
^He production in low mass stars often overpro- 
ducing ^He. Since this low mass behavior never- 
theless seems to be demanded by planetary neb- 
ula data, the data seems to favor models which 
have very little processing in general. 

We may also examine correlations between the 
calculated (pre) solar abundances of other pairs 
of light elements. Figure ^ shows D and ^He 
(pre)solar abundances. Here again, the largest 
segregation of points is due to the different initial 
D abundances for different r/. The shapes of the 
point distributions are similar for the different 
r] values, though the higher rj regions are more 
compressed. Note that the regions are shifted 
upwards for increasing rj, reflecting the higher Yp 
for these models. However, it is clear that the 
shift due to different Yp values is a much smaller 
effect than the spread due to chemical evolution 
effects, which give Ay anywhere from 0.04 to 
0.08. 

For a given r], the D-'^He correlations are 
not as immediately clear as for the D-^He plot. 
One can understand the trends by separating the 
points according to the '^He model features to 
be varied, namely the 8 — 12 Mq yields and the 
presence of mass loss in massive stars. For mod- 
els with no extra ^He yields (i.e. using just the 
adopted yields without mass loss and interpolat- 
ing between 8 — 10 Mq), the predictions appear 
as the lowest band of points; data for these mod- 
els alone appears in Figure |l^ (a). These models 
do indeed show an appreciable correlation, as we 
have argued they should: D depletion is the mea- 
sure par excellence of the amount of material pro- 

This is true for sufficiently early epochs and for sufficiently 
high gs, both of which are the case in our study. 
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cessed through stars, and '^He production is one 
consequence of this processing. Thus we expect 
and find a negative correlation for the "standard" 
case. 

Figure |l^ (b) shows the models with one 
or both of the helium enhancing features. We 
see here, like the D-^He plot, there are three 
"prongs" for each rj. These correspond (for the 
most part), to the three different cases of "non- 
standard" ^He yields. The lowest trend is for 
mass loss only, the intermediate (and some of the 
highest) prong is due to high 8 — 10 Mq yields, 
and the highest points come from a combination 
of the two. Here again, we see the negative corre- 
lation we expect. We also note that many models 
appear to fit the solar ^He. Indeed, while some 
are high and some low, a good number of both 
standard and nonstandard ^He yields give a suc- 
cessful fit. 

We should emphasize that the physics that un- 
derlies the D-^He correlation is different from 
that behind the D-^He correlation. In the lat- 
ter case, we have demanded that D be made into 
^He, and so we have explicitly built in the cor- 
relation between the two. The relation between 
D and ^He is somewhat less direct, as Galactic 
^He does not come from primordial D, but rather 
the destruction of D happens in processes which 
also produce '^He independently of the level of 
D. Finally, having found correlations between D 
and ^He as well as D and ^He, it one also ex- 
pects correlations between ^He and ^He produc- 
tion. These exist but are more complicated at 
there is an interplay of both the ^He and '^He 
evolutionary features. The basic trend, however, 
is unsurprising: a correlation between the two. 

Having examined the correlations between dif- 
ferent model calculations at the epoch of solar 
birth, we now turn to correlations at the present 
epoch. Here the underlying physics is the same 
and so we of course expect similar relationships 
between the model yields, but the observational 
constraints are different, and in the case of D, 
more severe. In Figure 15 we plot D versus ^He 



for the present epoch, i.e. we plot the abundances 
calculated for the ISM now. One notices the 
same features as the plot for the solar system, 
now smeared out with the intervening evolution. 
For ^He the observational situation is uncertain 
as the data shows a large dispersion; if one in- 
cluded the range of abundances spanned by the 
observations it would not constrain any of these 
points. This contrasts with the power of the 
presolar '^He, which offers a tighter constraint 
than the presolar D. 



The D observations in Figure 15 are from the 
HST observation of Linsky et al. (|1993 ), with 
the 2 — a variation shown. (Unfortunately, we 
are unaware of accurate ^He abundances in the 
ISM, and so we are unable to constrain our mod- 
els with this isotope.) If this number is to be 
taken seriously, it is a strong constraint on chem- 
ical evolution models. Here again, we find that 
the low-r/ points are very disfavored by the D ob- 
servations, with the lowest points in this regime 
needing to move at least a factor of 2 to meet 
the data. On the other hand, the spread in these 
points is larger than the presolar data, being al- 
most a factor of 3. Recalling the caveats above, 
we note that our models have trouble fitting the 
ISM D abundance given an initial D from a \ow-r] 
BBN model. 

The stringency in the observed presolar and 
ISM abundances for D, illustrated in figures 12 
and bears further investigation. As we have 
noted, D is the nuclide whose chemical evolution 
is the simplest, and thus we can expect a given 
model to calculate D evolution the most accu- 
rately of all the nuclides. It is thus interesting 
to ask what models survive the constraint that 
both the presolar and ISM D observations are fit 
to within 2 — cr. As is clear from figures 13 and 



15, the more stringent of the two cuts is comes 



from fitting the ISM observation. 

Demanding that models fit presolar and ISM 
D reduces the number of models from 1460 to 



231. In Figure 16 we plot the calculated presolar 
D versus ^He for these models. Comparing these 
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results to those of Figure we see that some 
of the models which could otherwise fit the (also 
stringent) presolar ^He are removed. However, 
many of these models still remain. Encourag- 
ingly, most of the models which successfully fit 
the presolar and ISM D also fit the (pre)solar 
^He and ^He. 

While it is not our purpose to find the "best 
models" for light element chemical evolution, it 
is useful to consider the reason for the success of 
the models which fit the observed (pre) solar D, 
^He, and ^He abundances, as well as the ISM D. 
Most importantly, these models all have a high 
rj and so begin with a very small amount of D 
and ^He— indeed, at D/Hp = 2.3 x 10-^ they 
are already smaller than the 2 — a upper limit 
to the presolar abundance. Thus these demand 
very little stellar processing of D to reduce it to 
the observed abundance, and therefore need to 
avoid producing much ^He. 

Interestingly, these models have very few other 
strong identifying characteristics. All low mass 
^He prescription are represented in an even way, 
as are all combinations of ^He prescriptions. The 
global characteristics include all SFR prescrip- 
tions and infall as well as closed box models. The 
two steepest IMF slopes are both represented as 
well, although the only mass limit range is the 
most restricted one (type 2 in table ^). 

We are left to conclude that D evolution is a 
strong constraint on chemical evolution models, 
but also that models meeting this constraint can 
still fit ^He and '^He. We have demonstrated that 
comprehensive evolution of the light elements is 
possible within our adopted framework of chemi- 
cal evolution, although this framework allows for 
substantial variation in its results according to 
the adopted model prescriptions and parameters. 

4.4. Additional Considerations 

44.1. D 

The program of observation of D in quasar 
absorption line systems, now just beginning, is 



likely to revolutionize how BBN is observation- 
ally constrained. The power and potential impor- 
tance and of this method impel a close scrutiny 
of its assumptions. In particular, it is implicitly 
assumed that the observed D abundance in high 
redshift systems is a direct record of the primor- 
dial abundance. However, by a redshift of z ~ 3, 
the universe it at least 1 Gyr old, and as much 
as 3 Gyr old (the range comes about since the 
t{z) relation depends on Oq)- Potentially, in this 
amount of time there might be some depletion of 
D in the protogalaxies that comprise the absorp- 
tion line systems; the question is how large such 
depletion could be. Also, it is noted that typical 
absorption line systems are not as metal poor as 
that studied by Songalia et al. (|1994| ), Carswell 
et al. ( |199^ ), and Tytler & Fan ( |T994| ); it is thus 
desirable to understand the evolution of D as a 
function of metallicity. 

Regarding the depletion of D in the absorption 
line systems, we note first that if these are proto- 
galaxies, we would expect them to take some time 
to form; however, we model the Galaxy once it 
has formed. Thus an object at a universal age of 
1 Gyr would correspond to a time earlier than 1 
Gyr on our models (indeed, Timmes, Lauroesch, 
& Truran ( 1995| ) argue for a significant delay be- 
fore the onset of galactic nucleosythesis). And 
furthermore, a glance at Figure |5| suggests that 
the D depletion at very early times is minimal. 

To address this question more systematically, 
we have considered the D depletion in all of our 
models for universes at a redshift of z = 3, where 
we have made the very conservative assumption 
that the protogalaxies form immediately and so 
time in our models corresponds to universal time. 
The correspondence between the observable z 
and the time t in our models depends on the 
adopted cosmology, in particular on the value of 
0. To be conservative we have chosen = 
(i.e. a curvature dominated universe, which has 
the largest age at z = 3), with a Hubble constant 
of Hq = 75 km/s/Mpc. For each of these cases 
we run our models to the time corresponding to 
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z = 3 and find the D abundance. To get a sense 
of how this correlates with the metal production, 
we plot the D abundance for these models against 
[Fe/H]. Results appear in Figure |l^ The Fig- 
ure confirms that to a high accuracy (> 20% and 
usually better than 5%) D may be considered un- 
depleted in these systems. Furthermore, Figure 
illustrates that for metal poor systems with 
[Fe/H] < — 1, and not necessarily at z = 3, the D 
depletion is minimal regardless of the IMF and 
infall choices. 

44.2. ^He 

The results from the previous section make 
clear that our ^He options can produce very dif- 
ferent ^He abundances at the birth of the solar 
system and in the ISM. We now turn to the issue 
of how these features affect ^He at low metal- 
licities, in particular how the model features af- 
fect the y— CNO relation studied in extragalactic 
H II regions. We show this relation for all four 
"^He model options (all at t] = 4x 10"^'^) in Figure 
p!g| . We see that the 8 — lOM© enhancement is ef- 
fective at increasing the ^He slope, making an en- 
hancement of about a factor of two over the stan- 
dard case. On the other hand, mass loss is unim- 
portant in this regime. This result is not surpris- 
ing when one considers that mass loss — at least 
as we have modeled it following Maeder (1992) — 
not only increases "^He but also CNO. Thus while 
the ^He production is enhanced, the tracers are 
as well and so the net effect on 1"— CNO is small. 
A note of caution, however, is that the mass loss 
calculations are difficult and so the results we 
use are subject to large uncertainty; also, they 
are only available for two metallicities. Improved 
modeling is crucial to help address the issue of 
the y-CNO relation. 

Figure (a) shows effect of secondary N evo- 
lution on the low metallicity 1"— N relation; this 
relation needs to be well understood to properly 
extrapolate from the observations to derive the 
primordial ^He abundance. Oftentimes a linear 
relation is assumed; we now estimate the error 



in assuming a linear Y— N relation for different 
cases of N evolution. For each model, we will try 
to reproduce the empirical fitting procedure by 
finding the slope of the 1"— N curve at a point and 
then extrapolating to get the estimated primor- 
dial ^He abundance Y^^^. We can then compare 
this to the actual Yp value in the model and thus 
compute the error 6Yp = Yp — 1^^* in this extrap- 
olation procedure. 

A potential complication to this procedure is 
that those models lacking features which enhance 
^He (mass loss and enhanced 8 — I2M0 produc- 
tion) will have smaller y— N slopes and so will 
be poor approximations to the extragalactic H II 
region data. Thus we will also keep track of the 
average slope one would assume for these regions 
if one did a linear fit; only the results with a 
sufficiently high slope will be admissible for this 
comparison. 

We perform the extrapolation at N/H = 50 x 
10"^ i.e. at a region surrounded by data and 
not at low enough metallicity to itself betray the 
low-N dropoff in ^He. The average slope is not 
the one used in the extrapolation but is taken at 
N/H = 100 X 10~^. This quantity serves only 
as a diagnostic to indicate the goodness of fit of 
the given model's y— N relation to the observed 
slope. As such, a central point better quantifies 
the average behavior over the whole region. 

In Figure |2^ we plot the error 6Yp in the 
^He linear extrapolation versus the average slope 
dY/dN one would assume for the region. Note 
the distribution in slopes; we will only consider 
models with slopes above 10 x 10^, which is a 
rough lower limit to the slopes derived from the 
H II region observations. For the models with 
large slopes, one can immediately discern the 
trends seen in Figure |^, here writ large. Namely, 
the models with no secondary N are well fit to a 
line, showing little error in the extrapolation. In- 
deed, the errors are fairly evenly spread around 
zero, with a width of about 0.001 units in the 
mass fraction. The models with some secondary 
and some primary N give 6Y ~ 0.002 — 0.004, 
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while models with purely secondary N can have 
corrections as large as 0.015. 

Thus the models with exclusively secondary N 
give 5Yp at the level claimed by Fuller, Boyd, & 
Kalen ( |T99lD , Mathews, Boyd, & Fuller ( |T99^ ), 
and Balbes, Boyd, & Mathews ( 1993| ). However, 
as pointed out by Olive, Steigman, & Walker 
( |1991| ) , the N-0 relation in extragalactic H II re- 
gions demands some degree of primary N; indeed. 
Figure fits the data well with a mixture of pri- 
mary and secondary N. For this case, the error 
in Yp extrapolation is small. Indeed, differences 
at this level are unimportant compared to the 
6Y ~ 0.008 observational uncertainty. 

It is amusing to note that linear extrapolations 
of H II region data often give slightly larger val- 
ues when using N as a tracer than when using O, 
with the difference being around 0.002 to 0.003 
units, just the level of the effect we see in the pri- 
mary -|- secondary N models. Clearly more work 
on this issue is needed, and the good agreement 
of our mixed N model may be somewhat fortu- 
itous, as we have not explicitly built models of 
the irregular galaxies in which these observations 
are made. Nevertheless, we our models suggest 
that the the effect of secondary N on ^He extrap- 
olation may not be as pronounced as originally 
suggested. 

5. Conclusions 

Our results on the sources of our model sen- 
sitivity confirm those of Tosi ( 1988| ), as we find 
that the elemental abundances of our models are 
most sensitive to the shape and mass range of 
the IMF. Indeed, the models we tried were in- 
compatible with a very shallow IMF slope for any 
mass range; this was the only model feature that 
was completely excluded by our solar system con- 
straints. Furthermore, we find that our "best" 
models, those fitting the solar CNOFe as well as 
solar and ISM observations of D, select an inter- 
mediate IMF slope, namely the Salpeter ( |1955| ) 
value. We also find significant sensitivity to the 



presence and timescale of infall. However, infall 
itself is strongly constrained by the low metallic- 
ity G-dwarf distribution. Finally, we do not find 
our results to depend strongly on the choice of 
the star formation rate. 

We also follow Tosi ( [1988| ) in finding that our 
results are very sensitive to the stellar nucleosyn- 
thesis yields employed. We emphasize in partic- 
ular the effect of systematic effects such as mass 
loss which are not included in many models. In 
our models mass loss can be an important source 
of ^He as well as CNO; more detailed calcula- 
tions of this affect are needed to address the im- 
portant question of the ^He contribution of high 
mass stars. We also find that our results are sen- 
sitive to the behavior of stars in the poorly un- 
derstood 8 — 12 M0range; this too provides an 
important uncertainty for ^He production. Work 
on these objects would be as useful as it will be 
difficult. Finally, we stress the need for a good 
understanding of the ^He yields in stars of all 
masses. While it would be quite useful to update 
the Iben & Truran ( 1978| ) calculation of ^He pro- 
duction in low mass stars, it is as important to 
understand ^He behavior in high mass stars, with 
emphasis on the effect of winds. 

Regarding the viability of light element chemi- 
cal evolution, it is an important point simply that 
we do find models which can fit solar and ISM 
constraints on D, '^He and ^He. This important 
"existence proof" demonstrates that the frame- 
work we have adopted, despite its roughness, is 
able to fit the data. Thus a consistent picture of 
light element evolution, from the early universe 
to the present day, can be drawn. 

From the perspective of BBN, it is reassuring 
that we find the effect of different ij so impor- 
tant for D and ^He chemical evolution. We see 
that the global chemical evolutionary models we 
have examined, while having a significant degree 
of model dependence, nevertheless preserve infor- 
mation about the primordial abundances. Thus 
we are encouraged in our use of solar and ISM 
data as meaningful constraints on BBN, with the 



26 



caveat that there remain significant uncertain- 
ties arising from chemical evolution. Indeed, as 
we have shown, for D and '^He the model depen- 
dences are larger than the observational uncer- 
tainties (for solar D and ^He as well as ISM D). 
Thus BBN constraints coming from the solar and 
ISM D and ^He are useful but should be regarded 
with caution. 

On the other hand, we do not find D to be 
subject to significant depletion or chemical evo- 
lutionary uncertainty at early epochs probed by 
measurements of quasar absorption line systems. 
We find the D abundances at this epoch to be a 
clean indication of the initial D abundance (and 
thus of rj) ; this result is independent of the partic- 
ular cosmology chosen. Thus more observations 
of these systems^] should eventually be able to 
put strong constraints on rj and thus on I^b- 

In contrast to the case of D and ^He, we find 
"^He evolution to be dominantly sensitive not to 
its initial abundance but to the uncertainties in 
its stellar production. However, we do find that 
within these uncertainties it is possible to find 
^He evolution at low metallicity which can fit the 
Ay/ AN and AY/ AO slopes better than stan- 
dard models would indicate. Further, we find 
that the low metallicity y— N trend is indeed sen- 
sitive to the prescence of secondary N. However, 
our models which best fit H II region N and O 
data require some degree of primary N and con- 
sequently do not lead to large deviations from 
linear y— N relations. This issue, however, mer- 
its further study. 

Future work in light element chemical evolu- 
tion in part would refine the framework adopted 
here, most importantly by the adoption of dif- 
ferent models of the halo phase. One might 
eventually hope to create a self-consistent stellar 
and chemical evolution scheme for the light ele- 
ments, along the lines of the approach of Timmes, 
Woosley, & Weaver ( |1995| ). However, there is ul- 



timately a need for a more fundamental under- 
standing of star formation and its relation to the 
dynamics and chemistry of the Galaxy. In par- 
ticular, an understanding of the physics behind 
the SFR and the IMF would be of tremendous 
value. Nevertheless, while uncertainties in chem- 
ical evolution remain, it is encouraging that they 
still allow models for the light elements to teach 
us about about BBN. 
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^Note however that deriving D abudances from the obser- 
vations may not be simple (Levshakov & Takahara 1995) 
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Table 1: Galactic Sources for the Light Elements 



Nuclide 


Production site 


Contribution to 






OBSERVED abundance 


D 


none 


not applicable 


^He 


low mass (< 3Mq) stars (?) 


> 50% (?)* 


^He 


stars of all masses 


~ 10 - 20% 



* ^He abundances in the ISM vary with mass of 



H II region 
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Table 2: Global Model Features 



iV±U LlCl i Cd If U.1 c 


41 

TT 




OX" XV 


1 

i. 






2 






Q 


2 






Ui CA.JJ I bill, 1 — / . O vjr J' 1 




O 


T — Cl-^rr 
T — LO vjryi 


liVir Dildpt? 


1 
± 






2 


from PDMF (Scalo 1986)* 




1 
± 


T — 1 n 




9 






o 


X = 1.7 


1 N/1 H vnocjc" IityitJ"C! 
IIVIX" llidoo lilllilo 


1 

± 






9 


(mi m ^ — (n ^ fin"! 




3 


imi m„) = fO 4 30) 


infall 




f = /oexp(-t/rinf) 




1 


/o = 




2 


Tinf = 2 Gyr; 99.9% atot 




3 


Tint = 4 Gyr; 50% atot 




4 


Tinf = 4 Gyr; 99.9% Utot 


N yields 




X^-^=Xl'{a + p{Xc/X^)) 




1 


a = 1,13 = 




2 


a = 0.5, /3= 1 




3 


a = 0, /3 = 1 




4 


a = 0,(3 = 2 


SN la normalization 


1 


A = 0.05 




2 


A = 0.007 



*Only calculable for exponential SFR 
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Table 3: Model Ranges for Solar D, ^He, and ^He 
Production 



V 




10^ D^in 


-Lu -L'max 


1f)^ iSfTp 


10^ •^Hp 


Y ■ 

mm 


Y 

^ max 


2 X 10" 


-iU 


5.90 


11.9 


4.26 


18.9 


0.250 


0.313 


4 X 10" 


-10 


1.82 


3.66 


2.00 


7.25 


0.258 


0.322 


6 X 10" 


-10 


0.912 


1.83 


1.44 


4.26 


0.262 


0.326 


V 






min 


(D+3He)^ax/(D+3He)^in 


Y —Y 

^ max ^ mm 


2 X 10- 


-10 


2.02 




3.03 






0.063 


4 X 10- 


-10 


2.01 




2.86 






0.064 


6 X 10- 


-10 


2.01 




2.58 






0.064 
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FIGURE CAPTIONS 

1. Big bang nucleosynthesis yields as a func- 
tion of the baryon-to-photon ratio r]. The 
*^He abudance is plotted as a mass frac- 
tion Yp] all other elements given as num- 
ber relative to hydrogen. Open points show 
abundance used for the three values of rj 
considered in this study. The dashed lines 
give the bounds for the Copi, Schramm & 
Turner (|1992 ) "sensible" bounds on rj. 



2. (a)The different adopted star formation rates, 
plotted as a function of time for models 
without infall. (b) The star formation rate 
ip = udtot plotted for the different models 
of infall. (c) As in (b), for -0 = i^crtot/^^- 

3. The (a) total and (b) gas mass surface den- 
sities clS cL function of time. Plotted for 
ip oc (Jtot with an IMF slope oi x = 1.35, 
with the four different infall prescriptions. 

4. Iron abundance function of time for 
different star formation rates. Data are 
from Edvardsson et al. ( |1993| ). 

5. For three different IMFs and SFR 5, evolu- 
tion of (a) D, (b) '^He and D+^He, and (c) 
^He. 

6. D and ^He evolution for the infall prescrip- 
tions. Note the enhancement of D and di- 
lution of ^He in the infall models. 

7. Low metallicities abundances of ^He, N, 
and O. data is for extragalactic H II re- 
gions, as reported by Pagel et al. ( |1992| ) 
and Skillman et al. (|99|). (a) N versus 
O. (b) The ^He mass fraction as a function 
of 0/H at low metallicities. (c) ^He as a 
function of N/H for different treatments of 
N evolution. 

8. (a) Plot of the yields for the full set of 
1184 candidate global models. Calculated 
for ?7 = 3 X lO"^'' and for a particular set 



of light element features. Abundances are 
calculated for the birth of the solar sys- 
tem and compared to observed solar abun- 
dances. Solar abundances are from Geiss 
( 1993| ) for D and ^He, and from Anders & 
Grevesse (1989) for all other elements; ±0.4 
dex limits to the solar abundances are in- 
dicated by the dashed lines, 
(b) Yields for the 267 candidate models 
having CNOFe within 0.4 dex of solar abun- 
dances. 



9. The G-dwarf distribution in two candidate 
global models. The model in (a) passes, the 
one in (b) does not; see discussion in text. 



Data is from Rana (1991). 



10. Solar abundances for the 2052 allowed global 
models with all light element model vari- 
ations. No cut has been made for CNO 
overabundances due to mass loss in mas- 
sive stars (see discussion in text). 



11. 



(a) As in figure 10, with the constraint that 
CNOFe be within 0.4 dex of solar; 1460 
models displayed. Note that imposing this 
constraint makes little change in the light 
element ranges. 

(b) As in figure |l^, for models with rj = 
4 X 10- 



-10 



12. The data of figure |TT| which reproduces 
CNOFe within 0.3 dex (short dashed line). 
For comparison, a level of 0.4 dex is shown 
in a long dashed line. 

13. Plot of D and ^He predicted solar abun- 
dances for the "final set" of 1460 models. 

14. (a) (Pre)solar D versus ^He for models in 
having no '^He enhacement features. 

(b) Solar D and ^He for models having ei- 
ther or both ^He enhacement features. 

15. Calculated abundances of D and ^He at the 
present epoch. As discussed in the text, it 
is unclear how to use the ISM ^He data as 
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a constraint; we have omitted it here. We 
show the D abundance in the ISM due to 
Linsky et al. ([1993D . 



16. Calculated D and ^He abundances at the 
solar birth. As in figure but models 
have been constrained to fit both the preso- 
lar and the ISM D abundances. 

17. D versus [Fe/H] for our models, plotted at 
a galactic age equal to universal time at 
z = 3. The redshift-time relation is derived 
for a cosmology with $7 = and Hq = 75 
km/s/Mpc. Dotted lines indicate primor- 
dial levels for the three adopted rj values. 

18. D versus [Fe/H] for the models of figure 
^ having different IMFs and for the high 
infall model 4 of table ^. Regardless of 
these model features, the D evolution at 
low metallicity is minimal. 

19. Mass fraction of '^He plotted as a function 
of (a) C, (b) N, and (c) O abundance in 
the low metallicity regime. Solid line: no 
^He enhancement; dashed line: 8 — 12 Mq 
^He enhancement; dot-dashed line: mass 
loss; dotted line, both enhancements. For 
r] = Ax 10"^°. 

20. The error 5Yp made in extrapolating the 
y— N relation linearly to determine the pri- 
mordial ^He abundance. Plotted for each 
model function of the average slope 
dY/dN of the Y— N relation at low metal- 
licity. Open circles are for models with only 
primary N, triangles are models with mixed 
primary and secondary N, squares are mod- 
els with low secondary N and hexagons are 
for models with high secondary N. Hori- 
zontal line indicates level of observational 
error; vertical line indicates approximate 
level of minimum slope for H II region data. 
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